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Abstract 

We theoretically study the Talbot effects resulted from the four- wave mixing and six- wave mixing signals, 
which are periodically modulated due to the coherence control effect. Corresponding to different dressing states, 
the enhancement and suppression conditions, that will affect the properties of the mult i- wave mixing signals are 
also discussed in detail. Such proposal can be useful in all optical-controlled pattern formation and propagation 
of light. 

1 Introduction 

Talbot effect, which was observed by H. F. Talbot in 1836 [1] and first analytically explained by Lord Rayleigh in 
1881 [2], is a near-field diffraction phenomenon that the light field spatially imprinted with periodic structure can 
have self-imaging at certain periodic imaging planes (the so-called Talbot planes). Such self-imaging effect holds a 
range of applications from image preprocessing and synthesis, photolithography, optical testing, optical metrology, 
spectrometry, to optical computing [3]. 

To date, studies on Talbot effect have been reported associated with atomic waves [4,5], surface waves [6,7], 
nonclassical light [8], or waveguide arrays [9]. Among the quiverful of researches, the nonlinear Talbot effect 
with second-harmonic (SH) were demonstrated both experimentally [10] and theoretically [11] in periodically poled 
LiTaOs (PPLT) crystal for the first time. The Talbot effect with beam imprinted with an electromagnetically 
induced grating (EIG) has been also proposed [12]. However, no progress in higher nonlinear optical process was 
reported before as far as to our knowledge. 

The crucial point and the guiding ideology to observe the Talbot effect is how to produce a spatial periodical 
incidence in transverse dimension. From this point of view, it is natural to realize that researchers have created 
spatially periodic four- wave mixing (FWM) and six- wave mixing (SWM) signals due to the periodic atomic coherence 
induced by standing waves (SWs) [13-16], which can play an important role in lasing without inversion [17], slow 
light generation [18], photon controlling and information storage [19-21], fiber lasers [22,23], and so on. The periodic 
pattern of the multi-wave mixing (MWM) signals can be flexibly controllable by adjusting the atomic coherence 
via changing the beam detunings. 

In this article, we investigate the MWM Talbot effect for the first time. And as the coexisting FWM and SWM 
processes have been observed in the Cd (S, Se) semiconductor doped glasses [24], the idea can be also executed 
in some solid crystals such as Pr-doped YSO crystals [24-27] besides atomic vapors. The article is organized as 
follows: in Sec. 2, we introduce the basic theory includes the energy level systems, some conceptions, and how to 
prepare the periodic FWM and MWM signals; in Sec. 3, we discuss the enhancement and suppression conditions, 
due to which the energy levels split; in Sec. 4, we investigate the Talbot effects of the periodic MWM signals; and 
in the last section, we conclude the article. 



2 Theoretical Model and Analysis 

We consider the FWM and SWM processes in the reverse Y-type atomic levels system derived by five beams, as 
shown in Figs. 1(a) and (b). A candidate for the systems are rubidium atomic vapors. Also, in Pr:YSO crystal, the 
energy-level system with |0) = | 3 # 4 (±3/2)), |1) = \ 1 D 2 (±3/2)}, |3) = \ 3 H A (±l/2)} can be used as a sub-system of 
our proposed system. Here we note that to guarantee the applicability of our proposal to all the possible candidate 



1 



systems, we choose parameters in simulation arbitrarily. Figs. 1(c) and (d) are the corresponding beam geometric 
configurations. The transition |0) — >> |1) with resonant frequency Q±q is probed by beam Ei (with frequency lji, 
frequency detuning Ai and wave- vector ki); |1) — > |2) with resonant frequency Q21 is pumped by beams E2 (with 
UJ2 1 A2 and k2) and E 2 (with co>2, A2 and k 2 ); |3) —> |1) with resonant frequency ^31 is pumped by beams E3 
(with A3 and k^) and E 3 (with CJ3, A3 and kg). The frequency detunings are defined as Ai = Q10 — lji, 
A2 = ^21 — ^2, and A3 = Qis — u)s • In the Cartesian coordinate introduced in Figs. 1(c) and (d), the wave- vectors 
are elaborated designed that Ei propagates ^-negative with an angle Q\\ E2 and E3 propagate along the opposite 
direction of ki , deviating from the ^-positive direction with 9\ ; E 2 and E 3 propagate symmetrically to E2 and E3 
with respect to z-negative direction. 

As shown in Fig. 1(a), if the fields E2, E 2 , E3 and E 3 are not strong sufficiently, we consider the undressed 
FWM signal E^ with ujp = ^1 and kp = ki + k3 — k 3 ; if E3 and E 3 are weak but E2 and E 2 are strong sufficiently 
to induce dressing effect (act as the dressing fields), we consider the singly-dressed FWM signal ~Efi with uj^i = oo\ 
and kpi = ki?. While if E3 and E 3 are also strong sufficiently to induce dressing effect, we consider the doubly- 
dressed FWM signal E^?2 with uop2 = oj\ and kp2 = kp. Even though E2 and E 2 will induce another FWM signal 
(P polarization), we still can get the expected FWM signal via a polaroid for its polarization selectivity. By using 
five beams with polarizations as shown in Fig. 1(b), a SWM signal E5 with ujs = ^1 and k# = ki +k2 — k 2 +k3 — k 3 
can be obtained. 




Figure 1: The schematic of the reverse Y-type atomic system to produce MWM signals, (a) kp represents the 
singly- or doubly-dressed FWM corresponds to a weak or strong k3, respectively, (b) Emission of the doubly-dressed 
SWM signal k#. The letters S and P in the brackets mean the polarization of the wave, (c) and (d) The Cartesian 
geometric configurations to generate FWM and SWM signals, respectively; The SWs along x and y are constructed 
by k2, k 2 and k3, k 3 , respectively. 

To obtain the MWM signal, we solve the density matrix equations and the Liouville pathways for p^ and p^ , 

the amplitude of which is proportional to that of the MWM signal. First, for the undressed FWM signal E^, we 

(3) 

can obtain p\ J via the Liouville pathway [28] as 

(0) Ei (1) (E 3 )* v (2) E' (3) 

Poo — * P10 > Pso — * Pio, (!) 

(3) -^FWM / x 

p ™ = ^p^' (2) 

where Gfwm = ^ioG3o(G 30 )* exp(zkp • r)exp(— iuit), d\ = Tio + zAi, ds = T3o + i (Ai — A3). The terms 
Gio — PmnEio/h and G' i0 = p mn E' i0 /h with i = 1,2 are the amplitudes of the Rabi frequencies of E$ and E^, 
respectively, where p mn is the dipole moment of transition between \m) and |n), which is driven by E^; and 
E' iQ are the amplitudes of E^ and E-, respectively. T pq represents the transverse relaxation rate between \p) and \q) 
for p 7^ and the longitudinal relaxation rate of \p) for p = q. It is obvious to see that the amplitude of p\ in Eq. 
(2) is spatial independent, so is uniform. 

Next, we consider the singly-dressed FWM process from E2 and E 2 , which are strong sufficiently. In the spatial 
interaction region, E2 and E 2 will interfere with each other and create a SW, which leads to periodic Rabi frequency 
amplitude \G2t(x)\ 2 = G\§ + G' 2 \ + 2G20C20 cos P (&2 sin#i) x). Therefore, the Liouville pathway for E^i will be 
modified into: 

(0) JSi (1) (E 3 )\ (2) K (3) /on 

POO > P\G 2t (x)±0) ^30 > P\G 2 t(x)±0)> W 

where |G2t (#)=!=) represent two dressed states produced by the spatial periodic dressing effect. According to Eq. 
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(3), we can obtain the spatial-dependent dressed density matrix element as: 



( x ) = ~ iGFWM (4) 

F1 [di + \G 2t {x)\>/d 2 ] 2 dz 

(3) 

where &2 = + i (Ai + A2). The amplitude of p F [ shows obvious periodic variation along x direction with a 
period of d x = A2/(2sin#i). 

When E3 and E 3 get strong significantly, we will get the doubly-dressed FWM signal E^. The Liouville 
pathway and the doubly-dressed density matrix element for this case are 

(0) Ei (1) (gsT (2) K (3) /-x 

POO > P\G 3t ±G 2t {x)±0) ^30 > P\G 3t ±G 2t (x)±0)i W 



= r , ^ " (6) 



,(3)^ _ — iGFWM 

" [di + |G 3 o| 2 M + |G2 t (z)| 2 M] 2 d 3 ' 

In Fig. 1(b), the spatially periodic dressing effect from E2 and E 2 remains, and the dressing fields E3 
and E 3 with same polarization (compared to the case of Eq. (6) with E3 and E 3 having orthogonal polariza- 
tions) can also interfere with each other and create periodic Rabi frequency amplitude \G3 t (y)\ 2 = G^ + G 3 2 + 
2G30G30 cos [2 (ks sin #2) y]. Therefore, this SWM will suffer from two dressing effects, which are both spatial pe- 
riodic with different periods d x and d y = A3 / (2 sin #2), respectively. The corresponding Liouville pathway and 
doubly-dressed density matrix element are 

(0) EV (1) (E 3 )* (2) K (3) (E 2 )\ (4) K (5) / ? x 

POO y P\G 3t (y)±G 2 t(x)±0) ~^ P30 y P\G 3t (y)±G 2t (x)±0) ~^ P2±0 > P\G 3t (y)±G 2t (x)±0) ' \ ( ) 

(5) _ iGswM / q \ 

Pl ° ~ [dl + |G 2 t(x)| 2 /d 2 + |G3t(^)| 2 /*] 3 d2d 3 ' 

where Gswm = G W G 20 (G f 20 )*G 30 (G f 30 )* exp(ik s • r) exp(-icjit). 

With these expressions, the spatial periodic pattern formation of MWM signals can be investigated. It is clear 
that the periodical variation of singly-dressed FWM, doubly-dressed FWM and doubly-dressed SWM signals with 
respect to x-axis at z = 0, are all derived from the periodical dressing effects. With different probe detunings, these 
MWM signals will show profiles with significant differences. So, in order to analyze such dependence, we consider 
the so-called enhancement and suppression conditions due to dressing effect before investigating the MWM Talbot 
effects. 



3 Suppression and Enhancement Conditions 

The characteristics of FWM are affected by the dressing effect, which depends on Ai, A2, and \G2t ( x )\ 2 • As shown 
in Fig. 2, because of |G?2t (^)| 2 7 |1) is split into two dressing states \G2t{x)±) with eigen-frequencies X\G 2t (x)±) = 
A 2 /2 ± v / A?,/4+ \G 2t (x)\ 2 . Therefore, \G 2t (x)±) are periodic along x. 

Absorption will be enhanced when the probe resonates with the dressing states, i.e., Ai = —\\G 2t (x)±)> which 
corresponds to the electromagnetically induced absorption (EIA) condition. Accordingly, the FWM signal will get 
enhanced resonantly. Thus we define Ai = —X\G 2t (x)±) as the enhancement condition. While when the probe 
reaches two-photon resonance (Ai + A2 =0), absorption will be suppressed (the EIT condition), and the FWM 
signal in such case will be suppressed correspondingly. Thus, we define Ai + A2 = as the suppression condition. 
So, if Ai is set at discrete values orderly from negative to positive and A2 is scanned, we can obtain the periodic 
enhancement/suppression condition and FWM signals along x. In Figs. 2(al)-(el), we display the periodic FWM 
signals along x corresponding to different Ai, and the insets are the dependence of the FWM signal intensity on 
A2 at z = and x = 0. In Figs. 2(a2)-(e2), we exhibit the corresponding split energy level diagrams, and in 
Figs. 2(a3)-(e3) the corresponding periodic properties of the split energy levels along x. In Fig. 2(al), Ai is 
relatively large and positive, and we only consider the enhancement because of the weakness of the suppression. 
When Ai decreases but is still positive, A2 will meet the enhancement condition (Ai = —\\G 2t ( x )+)) fi rs t, and 
then the suppression condition as shown in Fig. 2(bl). Because of Ai = in Fig. 2(cl), so that A2 can never 
meet the enhancement condition. The case in Fig. 2(dl) is opposite to that in Fig. 2(bl). And in Fig. 2(el), the 
enhancement is dominant again, the same as that in Fig. 2(al), for Ai is relatively large and negative. 
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Figure 2: The normalized intensities versus A2 and x of singly-dressed FWM with (al) enhancement effect, (bl) 
enhancement-suppression effect, (cl) suppression effect, (dl) suppression-enhancement effect, and (el) enhancement 
effect, corresponding to Ai = —5 MHz, —4 MHz, 0, 4 MHz, and 5 MHz, respectively. The insets are the incident 
intensity versus A2. (a2)-(e2) The corresponding split energy level (red lines marked |+) and | — ) on the right are 
split energy levels, the dotted line is the initial energy level) diagrams with dressing fields at x = 0. (a3)-(e3) The 
corresponding periodic split energy levels versus x. The other parameters are G 2 — 15 MHz, Tio = 5 MHz, and 
r 20 = 1 kHz. 



For the doubly-dressed MWM signals E^ 2 and E5, we take Ei? 2 to execute our discussions. First, the interference 
between E3 and E3 leads to periodic dressing effect, and therefore split the naked state into two dressing states 
\Gst(y)^=) with eigen-frequencies X\G 3t (y)±) — ~ A3/2 ± y/A|/4+ \Gst(y)\ 2 - E 2 and E 2 also act as dressing fields, 
so the first-order dressing state \Gst(y)^=) will split into second-order dressing states \Gst(y) ± G 2 t(x)=b) with 
frequencies 



-A 3 - VA| + 4|G 3t (y)| 2 , Aj L ±yg£+gK?ggP , Q . 

A \G 3t (y)-G 2t (x)±) = 7, 1 ' ^> 



-A 3 + VAl + 4|G 3t (t/)| 2 , A^ygf+ggggp 

A |G 3t (y)+G 2t (x)±> = 7, 1 7, ' ( W > 



where A' 2 = A 2 - {-A 3 - ^/A% + 4|G 3 *(y)| 2 }/2, A' 2 ' = A 2 - {-A 3 + ^/A% + 4\G 3t (y)\ 2 }/2, respectively. 

Similar to the case for E^i, we can also investigate the enhancement and suppression for E^ 2 and E5. And their 
enhancement and suppression correspond to Ai = —X\G 3 t(y)±G 2 t(x)±) an d Ai + A 2 = 0, respectively. The periodic 
doubly-dressed FWM signals along x corresponding to different conditions are shown by the mesh plots in Figs. 
3(al)-(il) with insets being the dependence of the signal intensity on A 2 at z = and x = 0.. The corresponding 
split energy level schematics are shown in Figs. 3(a2)-(i2), and the corresponding spatial periodic energy levels in 
Figs. 3(a3)-(i3). Doubly-dressed SWM enhancement and suppression effects are similar to the results shown in 
Fig. 3. In light of the SWM here is a two dimensional (2D) case, we exhibit two cases of the second order splitting 
energy levels as shown in Figs. 4(a) and (b). 

4 Talbot Effect of MWM Signals 

In the perspective of Fourier optics, the transfer function of a Fresnel diffraction system with z as the propagation 
axis, can be expressed as Hp(£) = exp (ik z z) exp (— inXz^ 2 ) in frequency domain [29], where £ is the spatial 
frequency and k z is the projection of k along z. The field of the MWM signal go (x,y) oc p(x,y), so it can be 
expanded into two-dimension (2D) Fourier series as go(x,y) = Em n=-oo Cm > n ex P J r m y/d y )) , and in 

frequency domain the equation above can be written as 

G {Z,v) = V ^sU-^-sL-^], (11) 



V d x J \ d 

nb,n=— 00 x 

where c m>n is the Fourier coefficient. So, considering the Fresnel diffraction, the MWM signal at a z distance is 

G (£, n) = Go (£, V) exp (ik z z) exp [-iw\ lZ (£ 2 + if)] ■ (12) 
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Figure 3: The normalized intensities versus A2 and x of (al) doubly-dressed FWM with enhancement effect, 
(bl) enhancement-suppression effect, (cl) suppression effect, (dl) suppression-enhancement effect, (el) suppression 
effect, (fl) enhancement-suppression effect, (gl) suppression effect, (hi) suppression-enhancement effect, and (il) 
enhancement effect, corresponding to A 1 = -18 MHz, -8 MHz, -5 MHz, -2 MHz, 0, 2 MHz, 5 MHz, 8 MHz, 
and 18 MHz, respectively. The insets are the incident intensity versus A2. (a2)-(i2) The corresponding split energy 
level diagrams with dressing fields at x = 0. (a3)-(i3) The corresponding periodic split energy levels versus x. The 
other parameters are A3 = 0, Gi = G3 = 15 MHz, Tio = 5 MHz, and = 1 kHz. 




Figure 4: (a) The surfaces represent \Gst (y) + Git (x) +) (top), |1) (grid), and \Gst (y) + Git (x) — ) (bottom), 
respectively, (b) The surfaces are \Gst (y) — Git (x) +) (top), |1) (grid), and \Gst (y) — Git (x) — ) (bottom), re- 
spectively. A2 = 2 MHz for (a) and —2 MHz for (b). The values of the other parameters are A3 = —3 MHz, 
G 20 = G' 20 = G30 = G' 30 = 15 MHz, Tio = 5 MHz, r 20 = r 30 = 1 kHz, Ai = 776 nm, A 2 = A 3 = 780 nm, 
0i =0 2 =0.3°. 
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Plugging Eq. (11) into Eq. (12), we end up with 



n 
d x 



Slrj 



m 



exp (ik z z) exp < — iirXiz 



(13) 



For simplicity, we let A2 = A3 and 0\ = # 2 , thus the periods along x and y are the same, i.e., d x = d y . Singly- 
and doubly-dressed FWM signals only concerns one SW formed by k 2 and k' 2 , and therefore we do not consider 
y component. In such a case, exp[— m\iz(n/ d x ) 2 ] = 1 if z = 2qd x /\\ (q = 1,2,3, •••), so after inverse Fourier 
transformation, we get 

9 (x) = go(x) exp (ik z z) . (14) 

Because of \g (x)\ 2 = \go(x)\ 2 , we can see the self-imaging of the MWM signals at z = 2qd x /\\, and zt = z\q=\ is 
the Talbot length. It is worth mentioning that the images on the zt/N planes are fractional Talbot images [30], 
where N is an integer bigger than 1. 




Figure 5: (al) The Talbot effect carpets for E^i under suppression conditions with Ai = 3 MHz and A 2 = 
—3 MHz. (bl) The Talbot effect carpets for ~Efi under enhancement conditions with A 2 = —3 MHz and — Ai = 
A 2 /2 + a/A^ + \G2t(x)\ 2 . (a2) and (b2) The intensity profiles at z = (black curves), z = z T /^ (red curves), 
z = zt/2 (blue curves), and z = zt (green curves), respectively. The other parameters are G20 = G f 2 o = 15 MHz, 
Tio = 5 MHz, and r 20 = r 30 = 1 kHz. 

We firstly choose E^i as the incidence to execute the simulation. As shown in Fig. 2, the spatial profile of 
the incident signal varies along A 2 with different Ai, and this variation leads to different diffraction process of 
FWM signal under different conditions. Figs. 5(al) and (bl) are the corresponding Talbot effect carpets under 
suppression conditions and enhancement conditions, respectively. Figs. 5(a2) and (b2) are the intensity profiles of 
the repeated images cutted at certain fractional Talbot lengths. From Fig. 5, we can not only clearly see the periodic 
singly-dressed FWM can reappear along z, but also see that the carpet stripes are oblique and such obliquity is 
more obvious under enhancement conditions. 

For Ei? 2 , the results are shown in Fig. 6. Clearly, we obtain the self-imaging of the incident E^ 2 at the Talbot 
plane again, no matter it is under the suppression or the enhancement conditions. In contrary to the case in Fig. 
5(al), which is obtained under suppression conditions, the Talbot effect shown in Fig. 6(al) seems oblique. Under 
enhancement conditions the carpet stripes are more obvious as shown in Fig. 6(bl), but the obliquity is almost the 
same with that shown in Fig. 6(al). By comparing the Talbot carpets shown in Fig. 6(bl) and Fig. 5(bl), we find 
the obliquity and the width of the stripes are almost unchanged. 

Last but not the least, we discuss the Talbot effect from E5. We consider two orthogonal SWs from two couple 
of dressing fields E 2 , E 2 and E3, E 3 simultaneously, to form a 2D lattice, which is periodic both along x and y as 
shown in Fig. 1(d), thus a 2D SWM signal will be excited. In Figs. 7 and 8, we first give the iso-surface plots 
of the Talbot effect of the 2D SWM signal under suppression conditions and enhancement conditions, respectively. 
And then we choose four intensity plots at certain places during propagation to show the details more clearly. We 
can see that at the Talbot length as shown in Figs. 7(d) and 8(d), the 2D SWM signals are reproduced. At half 
the Talbot length as shown in Figs. 7(c) and 8(c), the self-images shifted half period both along x and y. At 
one quarter of the Talbot length, fractional self-images can be seen as shown in Figs. 7(b) and 8(b), in which the 
images are twice as many as those in Figs. 7(a) and 8(a). Under enhancement conditions shown in Fig. 8, the 
images are more clear than those under suppression conditions as shown in Fig. 7. Here, we want to point out the 
reason that the periods of fractional Talbot effects shown in Figs. 5(b2), 6(a2) and (b2), and 8(b) seemingly not 
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Figure 6: (al) The Talbot effect carpets for Ei?2 under suppression conditions with Ai = 18 MHz, A2 = — 18 MHz, 
and A3 = 18 MHz. (bl) The Talbot effect carpets for ~Ef2 under enhancement conditions with A2 = —18 MHz, 
A3 = 18 MHz, and Ai based on Eq. (9). (a2) and (b2) The setup is as Figs. 5(a2) and (b2). The other parameters 
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Figure 7: The iso-surface plot is the Talbot effect for E5 under the suppression condition. The four panels are 
the contour plots of the Talbot effect at z = (a), z = zt/^ (b), z = zt/2 (c), and z = zt (d), respectively. The 
parameters are Ai = 18 MHz, A 2 = -18 MHz, A 3 = 18 MHz, G 2 o = G" 20 = G30 = G 30 = 15 MHz, Tio = 5 MHz, 
T20 = r 30 = 1 kHz, 



the half period at z = 0, is due to the inertial weakness of the numerical simulations. The Talbot effect changes 
dramatically during propagation, especially around small fractional Talbot lengths. Even a tiny deviation from the 
exact will bring seemingly imperfect quarter Talbot effect. Calculations with very high resolution will help us 
approach perfect fractional Talbot effects, but beyond our computer's memory. 
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Figure 8: The Talbot effect for E5 under the enhancement condition. The setup is as the Fig. 7, but with Ai 
according to Eq. (10). 



5 Conclusion 

In conclusion, we have studied the Talbot effects with completely controllable MWM signals. We have obtained 
the spatially periodic FWM and SWM signals by the interference between two dressing fields in a reverse Y-type 
atomic level system. The intensities of the MWM signals can be effectively controlled via enhancement conditions 
and suppression conditions derived from dressing effect. Talbot effect from singly- as well as doubly-dressed FWM 
is observed in the numerical experiment. Different from the case of FWM, for SWM, we construct a 2D lattice, from 
which the SWM is 2D modulated, to investigate the 2D Talbot effect. We find that the numerical simulations agree 
with the theoretical predictions very well. Our scheme is more advantageous in the controllability, compared to 
the previous studies, mainly attributed to the modulation of the MWM signal by the dressing effect. Our findings 
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not only enrich the understanding of the MWM theory, but also offer a different method to investigate the Talbot 
effect. 
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